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ABSTRACT 

The CoRoT observations of the {3 Cephei star HD 180642 uncover an unexpectedly rich 
frequency spectrum, in addition to several heat-driven modes. So far, two processes 
have been proposed to explain this behaviour: the presence of stochastic oscillations, 
and the excitation of time-dependent frequencies by nonlinear resonance. 

I argue for a third explanation for the observations, in the form of chaos due to 
the nonlinear behaviour of the dominant radial mode. 

The long-term frequency stability of the dominant radial mode is studied using 
archival data spanning roughly 20 years. Nonlinear time series analysis techniques are 
applied to the CoRoT observations, and the observations are compared with simula- 
tions of a simple nonlinear oscillator. 

I show that chaos offers one single explanation for many of the observed features, 
such as the structure in the autocorrelation of the power spectrum, the long-term 
frequency shift, the power excess and the wide range of frequencies in the power 
spectrum. However, the mixture of opacity-driven linear oscillations and nonlinear 
oscillations complicate the nonlinear time series analysis techniques. 
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1 INTRODUCTION 

HD 180642 is a large-amplitude (3 Cep star observed with the 
CoRoT space satellite Fridlund et al. (j2006j). Its rich fre- 



quency spectrum ( Degroote et al.||2009 ) immediately led to 



, debate on the interpretation. Belkacem et al. (2009) inter- 
preted the power excess at high frequencies (130 — 300 /xHz) 
as due to solar- like oscillations. [Degroote et al. (2009), how- 
ever, noted that there was also power excess at lower fre- 
quencies, and focused on the nonlinear resonances of the few 
high-amplitude modes and their effect on the high-frequency 
regime. These were subsequently confronted with evolution- 
ary models by Aerts et al. (2011). In the latter paper, the 



left-over frequencies were interpreted as being excited by 
nonlinear resonance with the dominant, high- amplitude ra- 
dial mode. 

Aside from the high-quality CoRoT observations, 



archival ground and space based data (Aerts 2000| Uyt- 



terhoeven 2008) span roughly 20 years, allowing to study 
long term temporal stability of the detected modes. Pe- 
riod changes are observed in various types of pulsating 
stars, in various stages of their evolutionary stages. Such 



changes were discovered in white dwarfs (Winget et al. 



1991 ) and (5 Set stars ( 


Breger & Pamyatnykh||1998 


). Pigul- 


ski & Pojmanski (2008 


) observed a frequency decrease {D ^ 



ASAS data, suggesting that period changes might occur in 
other /3 Cep stars as well. The large scatter in the 0-C dia- 
grams of Cepheids ( Szabo et al.|2011 ) hints for variability in 



the dominant frequency, while Kolenberg et al. (2010) finds 



smooth modulations of the main frequency in RR Lyr stars. 
Finally, Zijlstra & Bedding (2004) identifies different types of 
period changes, from sudden changes or irregular patterns, 
to smooth modulations in Mira stars. Bar the two first ex- 
amples, most of the observed period changes are larger than 
expected purely from evolutionary effects, or even have the 
wrong sign. This can possibly be attributed to 'evolutionary 
weather' (D. Kurtz, priv. comm.), in which case the evolu- 
tion of the star on the main sequence is not simply described 
by a continuous expansion of the star. Breger Pamyatnykh] 
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( 1998 ) propose two additional origins of period changes: the 
existence of a stellar companion and nonlinear mode inter- 
actions. The latter explanation is definitely appealing in the 
case of large amplitude pulsators. Nonlinearity can also give 
rise to period doubling events (as is observed in RRLyrae 
star, Kolenberg et al. j2010), and chaotic behaviour. The 
latter comes in many different forms, from erratic variabil- 
ity without a clear periodic signature, to mild forms of chaos 
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where a period signature is clearly detected, but not every 
cycle is exactly reproduced from one to the next. This can 
introduce a characteristic signature in the Fourier spectrum. 

In this paper, I report on the detection of a period 
change in the /3Cep HD 180642, for which over 20 years of 
photometric observations are available, among which a 150 d 
continuous, high quality CoRoT light curve. Aside from the 
apparent period change, also power excess is observed. As 
third alternative for the interpretations given by |Belka-| 



cem et al. (2009); Aerts et al. (2011), I show that chaotic 



behaviour can explain both the period change and power ex- 
cess. In Section 2, I present the archival photometric data, 
which is used in Section 3 to derive the magnitude of the 
frequency shift of the dominant radial mode. Subsequently, 
the power excess is studied and compared with simulations 
of a simple nonlinear oscillator. Nonlinear time series anal- 
ysis techniques are then applied in an attempt to firmly 
establish the presence of chaos. 
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Figure 1. Overview of photometric time series of HD 180642 
(blue denotes Hipparcos data, yellow is data from the Swiss 70 cm 
telescope, green from Mercator/P7, red from AS AS, cyan from 
Konkoly, magenta from SPMSNO and black from CoRoT. 



Table 1. Archival photometric datasets used in this paper. 



2 OBSERVATIONS 

The Hipparcos sateUite observed HD 180642 from 12 March 
1990 for almost three years. Since then, 191 observations as- 
sembled with the photomultiplier P7 attached to the 0.7m 
Swiss telescope at La Silla and to the 1.2m Mercator tele- 
scope at La Palma were made and added to the dataset. 
Next, the data from two observing campaigns at the Konkoly 
Observatory in Hungary and at the National Observatory in 
Sierra San Pedro Martir in Mexico were added. Moreover, 



there are 379 archival ASAS data points (Pigulski&Pa 
jmahski||2008|). F inally, I add the CoRoT light curve (see 
Auvergne erar]|2009l [Degroote et al.||2009[ for more de- 
tails) , containing more than 300 000 single observations over 
a time span of about 150 d, with far superior precision than 
the previous datasets. This brings the total time span of 
observations to almost 20 years (Fig.[T]). The characteristics 
of the different datasets (year of first measurement to, to- 
tal time span T and number of observations n), as well as 
the first frequency value / determined for each of these data 
sets are summarized in Table[l]and shown in Fig. [2] (top left 
panel). The frequency from the Konkoly observatory devi- 
ates strongly from the others, nearly simultaneous data sets. 
The reason is the awkward window function of this dataset, 
which leads to an alias of the true frequency having the 
highest amplitude in the power spectrum. 



3 DATA ANALYSIS 
3.1 Period change 

The evolution of the dominant period has already been ex- 
plored briefly by (Degroote et al. 2009), but their analy- 



sis was not quantitative and no physical interpretation was 
given. 

To quantify the decrease, I assume a simple linear de- 
crease and apply an adapted phase dispersion minimization 
(PDM) method (Stelhngwerf 1978] |Cuypers| 



1986) to the 



combined datasets, omitting the CoRoT dataset. The lat- 
ter is only used to check the findings, since the high cadence 
and relatively short time completely dominate the frequency 
search procedures. I test for 1000 trial frequency shifts D 
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between -4.5 x 10~^ d~^ and 0.5 x 10~^ d~^ in steps of 
5 X 10~^^ d~^ and calculate at what frequency the PDM test 
statistic (G) reaches its minimum (Fig. [5] top right panel). A 
clear minimum is found around Di — (—2.3=1=0.2) x 10~^ d~^ 
(a period increase of 2.4 seconds per century), which corre- 
sponds to a reduced scatter in the phase diagram compared 
to Di = Od"^ (Fig.[2] bottom panels). The error on the pa- 
rameter is determined via the covariance matrix obtained 
from the nonlinear optimization routine, and are confirmed 
with those obtained from an F-test applied on the statis- 
tic in an interval around the parameters. 

Excluding the Hipparcos measurements decreases the 
time span of the observations but yields consistent results 
(L>2 = (-3.0 ±0.5) X 10"^ d"^). Including the entire CoRoT 
dataset also yields consistent results within a 3cr confidence 
interval {D^ = (-2.61 ± 0.06) x 10"^ d"^). Regardless of 
the included datasets, there is an increase in variance re- 
duction and a decrease in AIC and BIG criterion (for a def- 
inition, see, e.g., Degroote et al. 2009|) when allowing for 
a frequency shift, confirming that a linearly decreasing fre- 
quency gives a better description of the data. The sparseness 
of the data sets, however, make it difficult to assess if the 
frequency shift is linear, quadratic, periodic or following yet 
a different pattern. The simple conclusion is that assuming a 
variable frequency, the data is significantly better described 
than without it. 

Possible causes for artificial frequency shifts between 
two datasets well separated in time, can be different time 
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conventions (e.g, Aerts, Christensen-Dalsgaard, & Kurtz 



2010 , Chapter 4). However, the results are insensitive to 
smah errors in the calculations of the time points. Deviations 
as large as 10 s translate to a shift in the phase diagram of 
only 0.06%, which has no influence on the results. Also, the 
results obtained from the different datasets from Hipparcos, 
CoRoT and ASAS separately, show a frequency decrease 
which is independent from any mismatch between different 
time coordinate systems (e.g., Julian Date or Modified Ju- 
lian date). A two-sample t-test shows that the probability 
that the frequencies in the Hipparcos and ASAS data are 
equal is below 0.1%. Finally, the aliasing problem which is 
present in the Konkoly data set, is not present in the Hip- 
parcos and ASAS data sets, eliminating the possibility that 
aliasing problems are at the basis of the frequency differ- 
ences. Thus, the observed frequency shift must be intrinsic 
to the target. 

A first physical explanation for period changes is light- 
time effects due to binarity. Light-time effects predict a ra- 
dial velocity of 



translating to a radial velocity change of ^ lOkms"^ over 
the 20 years of photometric observations of HD 180642. The 
datasets currently available are not sufficient to reveal pos- 
sible modulations. However, there are no signs of binarity in 
the high signal-to-noise, long term spectroscopic campaign 
done by Briquet et al. (2009), putting tight constraints on 



the luminosity of a hypothetical secondary object. 

I confront the period increase found in the combined 
datasets with expected values from a grid of evolutionary 
models (see also Aerts et al.|2011 ) calculated with the CLES 
( Scuflaire et al.|2008 ) code, where the adiabatic approxima- 
tion is used to calculate the frequencies with the LOSC code 
(Scuflaire et al. 2008) . The evolutionary frequency shift was 
calculated by extracting all ^ = modes from the grid, and 
taking the numerical 3-point Lagrangian derivative along 
the evolutionary track. As the star slowly expands during 
the main sequence, its mean density decreases together with 
the frequencies of the radial modes. The observed frequency 
decrease is at least 10 times larger than models predict, ex- 
cept for stars in the hydrogen shell burning phase. Previous 
seismic modelling excludes these stages, and given the short 
time scale of the hydrogen shell burning stage, it is unlikely 
to observe a star during this stage. 

After arguing against an instrumental, evolutionary and 
external explanation for the frequency shift, I propose an- 
other possibility: nonlinear effects. Although nonlinear ef- 
fects are not typically expected to produce consistent long 
term frequency variations, slight changes in phase timing can 
mimic long term trends when the observations are divided 
in distinct groups containing large gaps. In the following sec- 
tions, I look for more signs of nonlinear ity, and argue that 
this hypothesis can provide an alternative explanation also 
for the power excess observed in the frequency spectrum. 



3.2 Power excess 

A standard frequency analysis of the CoRoT light curve 
with iterative prewhitening results in thousands of formally 
significant peaks. Such a model is not very meaningful in 
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Figure 2. Upper left: Frequency determinations of ground and 
space based observations (same colour coding as in Fig.Q. Ver- 
tical bars denote la error in frequency, corrected for correlation. 
Horizontal bars denote time span for the frequency determina- 
tion, circles denote midpoint of observations. The offset on the 
frequency from the Konkoly data is a result of a bad window func- 
tion. Upper right: minimum PDM test statistic as a function of 
trial frequency shift. Lower left: phase diagram with D = Od~^. 
Lower right: phase diagram with D — 2.6 x 10~^ d~^. 



view of the limited frequency resolution. In contrast to 
the fairly constant-amplitude background excess in the fre- 
quency spectrum observed in HD 50844 ( [Poretti et al.|2009 ), 
the spectrum of HD 180642 is highly structured, even af- 
ter prewhitening 500 frequencies. [Kallinger &; Matthews] 
( 2010) suggested granulation to be the cause for the back- 
ground power observed in two A-type stars. In the case of 
HD 180642, that model does not give a satisfactory fit to 
the structures. 



Smolec &: Moskalik (2007) suggested that 



amplitude saturation could play an important role in excit- 
ing hundreds of resonant modes in /3 Cephei stars, perhaps 
even into the gravity mode regime. The number of observed 
frequencies in HD 180642 is, however, much larger, and the 
range of frequencies is large compared to the range predicted 
by excitation calculations. On the other hand, ( |Buchler"et 



al. 1996) interpreted the multiperiodicity in the RV Tauri 



star RScuti as low-dimensional chaos, and not as actually 
excited modes. 

In the following, I investigate if the nonlinear behaviour 
of the dominant radial mode could play a role in the power 
excess in the frequency spectrum. Confirmation of such non- 
linearity could also possibly explain the observed shift in 
frequency as due to slight irregularities in the dominant fre- 
quency: if the dominant mode is suffering from chaos, it is 
expected that not every pulsation cycle repeats itself exactly, 
which can mimick long-term frequency or phase shifts in 
grouped observations such as the archival data set at hand. 

It is known that even simple oscillators, e.g. those that 
are driven and damped, can show non-repetitive behaviour, 
as long as the amplitude of oscillator is in a particular range. 
Tanaka & Takeuti ( 1988 ) proposed and investigated a sim- 



ple ID model for nonlinear pulsations, which was used for 
simulations in Votruba et al. (2009). There, it was shown 



that chaotic oscillations can introduce spurious periods in 
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the periodogram. Their nonadiabatic oscillator is described 
by 



d^x{t ) jix{t) 



dt 



- uj x{t) — z(t), 



(1) 



where x{t) describes the surface displacement, /ll is the driv- 
ing factor, Lu is the frequency and z(t) the driving force. 
iTanaka & Takeutil (|1988|) set 



dz{t) 
dt 



dx(t) 
dt 



-pz{t) 



dt 



dt 



(2) 



The parameter p mimics the energy lost during a pulsation 
cycle, e.g., via radiation, while the nonlinear term preceded 
by the parameter s introduces a phase shift between the 
cooling and heating. The system can be solved using a a 
standard iterative Runge-Kutta method. As an example I 
set — 0.5, II — —0.5 (driving), p — 4.0 (radiative cool- 
ing), and s = 1.0 (phase shift between cooling and heating) 
by default, and only change how the material responds to 
volume changes via the parameter f3. High {3 values result 
in a strong drop of driving force when the material expands, 
in which case the amplitude of the oscillation should be low 
and thus the linear approximation is valid. Indeed, a value 
oi 13 — 2 gives a monoperiodic oscillator with a sinusoidal 
phase shape. On the other hand, a value of /3 = 0.7 drives the 
amplitudes to much higher values into the nonlinear regime, 
resulting in a distorted phase shape. 

To mimic the behaviour of HD 180642, 1 start from a set 



of parameters used in Tanaka & Takeuti (1988): a — A^h — 



^^IJi — — 0.5,0^ = —0.7, and s = 5. Changing the value of 
p between —0.04 and —0.08 creates a Feigenbaum diagram, 
that allows to carefully pick a value that restricts the size of 
the attractor in phase space, as is observed for HD 180642. I 
choose a value p — —0.726, emphasizing that the setup of the 
simulation is only chosen to qualitatively study the influence 
of nonlinear behaviour on the frequency spectrum. The aim 
is not to quantitatively model the dominant frequency of 
HD 180642 using Eq.[2] 

I generate a time series of 175 200 points, add Gaussian 
noise with a standard deviation similar to the CoRoT data, 
and apply standard iterative prewhitening tools as it was 
done for the stellar data in ( Degroote et al.|2009 ) to the sim- 
ulations. In Fig. [3] the original frequency spectrum is shown 
(top panel), where the harmonics and half- frequency (from 
the period-doubling phenomenon) are clearly visible. Some 
power excess is visible over the entire range of harmonics. 
Iteratively removing 13 frequencies makes the power excess 
at low amplitudes better visible. In a similar way as was 
found for the CoRoT light curve of HD 180642, the num- 
ber of frequencies from an iterative prewhitening analysis 
amounts to several hundreds to thousands. In contrast to 
red noise components or granulation signal, the background 
amplitude spectrum does not appear as a plateau, but rather 
as a group of dense peaks. 

The previous exploration shows that the frequency 
spectrum is contaminated by a huge number of high signal- 
to-noise ratio peaks, which are not connected to eigenfre- 
quencies of the system (there is only one in the simulated 
system). An autocorrelation (Fig.[4| of the frequency spec- 
trum reveals structure, which is similar to what is expected 
from stochastic pulsations (e.g., [Carrier Bourban|[2003t, 
and remarkably similar to the one shown in |Belkacem et al.| 
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Figure 3. Original frequency spectrum of a simulated nonlinear 
oscillator containing chaos (top panel), after removing 13 frequen- 
cies (middle panel) and after removing 300 frequencies. 
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Figure 4. Autocorrelation of residual frequency spectrum. 



( [2009 [, see the supplementary material). In the next Section, 
I attempt at quantifying to what extent HD 180642 exhibits 
nonlinear behaviour. 



3.3 Nonlinear analysis 

Chaotic signals can be recognised by their characteristic 
phase space behaviour. While monoperiodic oscillators have 
only one orbit in phase-space (i.e., each cycle repeats ex- 
actly), chaotic attractors appear in the case of the chaotic 
oscillators. 

When the dynamics of the system is known, the phase- 
space construction is straightforward and chaotic behaviour 
can be readily determined. In the case of observed time 
series, however, a different approach has to be followed. 
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Deterministic chaotic signals can be characterised by the 
largest Lyaponuv exponent (Ai), which is the average ex- 
ponential rate of divergence (or convergence) between two 
nearby points in phase space ( Wolf et al.|1985' ), 

d(t) = Cexp(Ait), 

where C is a normalisation constant. Lyapunov exponents 
thus deliver a measure of how quickly two neighbouring 
points in phase space separate. Strictly periodic signals have 
Ai = 0, because the orbits of any two points are the same. 
When Ai < 0, the orbits converge (e.g. in the case of a 
damped oscillator), and if Ai > 0, the orbits are unstable 
and chaotic. The basic steps to determine the maximum 
Lyapunov exponent thus consist of reconstructing phase- 
space and measuring the rate of divergence. I follow the 
false neighbour method approach of 'Rosens tein, Collins, fc| 
|de Luca (1993) , which is fast, easy to implement, and reliable 
also in the case of small datasets. Extensive numerical tests 
on systems of different chaotic degrees (i.e systems for which 
different embedding dimensions are required) as well as light 
curves with a highly multiperiodic signal, have shown that 
also the latter require high embedding dimensions, render- 
ing nonlinear time series analysis techniques on light curves 
containing highly multiperiodic signal less useful. In Fig.|5] 
the influence of the choice of embedding dimension is shown 
on the number of false neighbours, as is the influence of 
noise. For monoperiodic, the embedding dimension is very 
low, since the number of false neighbours drops quickly to 
zero. In highly nonlinear systems, the number of false neigh- 
bours drops only significantly below 100% at high values of 
the embedding dimensions. However, we see the same be- 
haviour for highly multiperiodic (but linear) systems. 

I conclude that, although it is very promising to sepa- 
rate monoperiodic signals and noise from chaotic time series, 
this is not the case for highly multiperiodic signals. Based 
only on the method of false neighbours and in the presence 
of noise, the same behaviour is found for the Lorenz II equa- 



100 



tions ( Lorenz|1984[ Kennel, Brown, &: Abarbanel|1992 ) and 
the multiperiodic signal. The limiting factor for the separa- 
bility of multiperiodicity and chaos will thus be the noise 
level and frequency resolution (or the number of cycles ob- 
served) . 

Choosing embedding dimensions between de — 2 and 
de = 7 and an embedding delay of r ^ 0.02 d (the value 
where the autocorrelation function drops below 1 — 1/e, 
Rosenstein, Collins, de Luca||1993 ), there is a high expo- 



nential separation rate, which reaches saturation level dur- 
ing the short time span of one orbit. This result is not sen- 
sitive to the number of modes we prewhiten beforehand to 
minimize the influence of multiperiodicity, nor to the choice 
of the embedding delay r. I conclude that the frequency 
spectrum of HD 180642 is dominated by the signal from ac- 
tual pulsation modes, so that nonlinear time series analysis 
techniques have difficulties detecting chaotic behaviour un- 
ambiguously in the CoRoT light curve of HD 180642. 



4 CONCLUSION 

In this paper, I proposed the presence of nonlinear behaviour 
as an alternative explanation for much of the low-amplitude 
power excess at low and high frequencies in the CoRoT light 
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Figure 5. Determination of embedding dimension using the false 
neighbour method (*=nonlinear nonadiabatic oscillator (NNO) 
with /5 = 2, diamond: NNO with /3 = 0.7, x =NNO with ^ = 1 
plus: NNO with /3 = 0.8, bu llet: Lorenz II ( [Lorenz 1 1984] [Kennel,] 
[Brown, Abarbanel[[T992| , black square: multiperiodic signal, 
black triangle: pure noise). Since multiperiodic signals also require 
a high embedding dimension, the false neighbour method cannot 
distinguish these types of signals from true chaotic signals. 



curve of HD 180642. The first explanation by [Belkacem et[ 
al. ( 2009 ) invoked the presence of solar- like oscillations, but 



failed to explain the lower frequency power excess. The sec- 



ond explanation was given by Aerts et al. (2011 ) via nonlin- 



ear resonance, but requires the presence of a dense spectrum 
of eigenmodes, containing many thousands of frequencies. 

It was not possible to unambiguously determine the ori- 
gin of the power excess, but from the point of view of model 
simplicity, I argue in favour of chaos, since it explains much 
of the observed behaviour without invoking new physical 
processes. At least four anomalies in the observations of 
HD 180642, compared to the expected relative scarce and 
structured frequency spectrum delivered by linear theory, 
are explained by the chaos hypothesis. First, it explains 
the observed period change, which is not compatible with 
current stellar evolution models, as a natural side effect of 
nonlinear behaviour. This is perhaps also applicable to the 
observed smooth period changes observed in RR Lyrae stars 
by^Kolenberg et al.j (2010J or the scatter observed by ^Szabo] 



et al. (2011 ) in Cepheids, or by Zijlstra & Bedding (2004) in 



Mira variables. All of these variables show large amplitude 
variations, driven far into the nonlinear regime. Second, it 
explains at least qualitatively the signal in the autocorre- 
lation. In contrast to red noise components or granulation, 
chaos can introduce structure in the form of dense forests of 
peaks in the power spectrum, introducing spurious peaks in 
the autocorrelation diagram. Third, it explains the wealth 
of peaks in a wide frequency range, that are not necessarily 
eigenfrequencies. Finally, if all the frequencies in the Fourier 
spectrum are true frequencies excited in the star by any 
mechanism, the amplitude saturation calculations of [Smolec| 



& Moskalik ( 2007 ) would most likely predict a lower satura- 



tion amplitude in the case of HD 180642. Their results show 
that the saturation amplitude roughly scales as n~°"^ with 
n the number of excited frequencies, and decreases even fur- 
ther with lower metal abundances. The observed visual am- 
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plitude of HD 180642 is at the high end of their calculations 
with collective saturation (their Fig. 7), while the number 
of observed frequencies in the pressure and gravity mode 
regime is at least an order of magnitude larger and the ob- 



served metal abundance is significantly sub-solar ( [Briquet et 
aL] [2009t . The presence of artificial frequencies in the light 



curve arising from the chaos, would strongly decrease the 
true number of modes taking part in the saturation effect, 
making the CoRoT observations more in line with the pre- 
dictions of ISmolec & Moskalikl (|2007|). 



A possible way forward from the observational point 
of view would be to continue long term monitoring of the 
dominant mode of HD 180642, which would allow us to check 
if the evolution of the frequency follows a consistent trend 
or is rather of chaotic nature. Most likely, the GAIA satel- 
lite will observe this star at least occasionally for an ex- 
tended period of time, which will, dependent on the launch 
date, add at least another ten years to the time span. The 
GAIA satellite will then also observe other high-amplitude 
(3 Cep stars, which would allow us to check if these period 
changes are common in /3 Cephei stars, even for low ampli- 
tude modes. If no correlation with amplitude can be found, 
the idea of short-term evolutionary fluctuations ('evolution- 
ary weather') should be favoured. From theory, it could be 
worthwile to apply specifically tuned nonlinear hydrody- 



namical models as presented by Smolec & Moskalik (20071, 



to reproduce the observed amplitude of the dominant mode. 

Still, the nonlinear timeseries analysis showed that the 
nonlinear behaviour, if present, is definitely not strong in 
the CoRoT data. The amplitude of the power excess is much 
lower than the amplitudes of the strongest modes (especially 
the dominant radial mode), and most of the high signal- 
to-noise peaks are due to additional pulsation modes. This 
is also confirmed by the successful modelling by [Aerts et] 
al. (2011|. We cannot, however, confidently identify which 
peaks of the low amplitude signal is due to pulsations, and 
which are more likely to arise from the nonlinear behaviour. 
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